Molecular Tilt on Monolayer-Protected Nanoparticles 
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The structure of the tilted phase of monolayer-protected nanoparticles is investigated by means of 
a simple Ginzburg-Landau model. The theory contains two dimensionless parameters representing 
the preferential tilt angle and the ratio e between the energy cost due to spatial variations in the 
tilt of the coating molecules and that of the van der Waals interactions which favors uniform tilt. 
We analyze the model for both spherical and octahedral particles. On spherical particles, we find 
a transition from a tilted phase, at small e, to a phase where the molecules spontaneously align 
along the surface normal and tilt disappears. Octahedral particles have an additional phase at 
small e characterized by the presence of six topological defects. These defective configurations 
provide preferred sites for the chemical functionalization of monolayer-protected nanoparticles via 
place-exchange reactions and their consequent linking to form molecules and bulk materials. 



Introduction. — Monolayer-protected nanoparticles 
(MPNPs) are small aggregates of metallically-bonded 
atoms 2-8 nanometers in diameter, coated with organic 
ligands. Organic materials readily absorb on surfaces of 
metal and metal oxides to lower the interfacial energy be- 
tween the substrate and the environment. These adsor- 
bates naturally alter the surface properties and can have 
a significant influence on the stability of nanostructures 
[U [2]. MPNPs have received considerable attention in 
the past few years due to the unique chemical and physi- 
cal properties that make them promising candidates in a 
wide array of potential applications, including electronics 
[3l , nanomedicine [5H7j , optics ^ , chemical microfab- 
rication [9 and biology p!Q] . 

Gold nanoparticles coated with alkanethiol chains, e.g. 
CH3(CH2)mSH, form one of the most thoroughly investi- 
gated classes of MPNPs because of their ease of fabrica- 
tion. Due to the strong interaction between alkanethiols 
and the gold substrate, the sulfur head groups of thi- 
ols chemisorb on specific sites forming a commensurate 
lattice on the substrate. At low surface coverage the two 
free lone pair electrons of sulfur drive the alkyl chain close 
to the gold surface (lying down phase). As the coverage 
increases van der Waals (vdW) interactions between the 
alkyl chains become important and the resultant compe- 
tition leads to a configuration with the molecules stand- 
ing upright on the substrate, but tilted with respect to 
the substrate normal. The minimum energy packing is 
achieved with the chains tilted at specific angles dictated 
by interlocking of the zigzag skeletal structure of the car- 
bon atoms (Fig. [T]). 

The morphology of a monolayer of tilted alkanethiols 
on the faceted surface of a gold nanoparticle has been 
the subject of intense investigation for over ten years 
pTHT4] . By means of molecular dynamics simulations 
it was found that, as a consequence of the seamlessness 
of the underlying substrate, the alkane chains tend to tilt 
uniformly across the nanoparticle, rather than assuming 



their equilibrium tilt angle on each facet, which would 
generate an energetically costly mismatch at the edges. 
On a closed surface such as a sphere or a polyhedron, 
however, perfectly uniform tilt can never be achieved. 
The projections of the tilted alkane chains on the tangent 
plane of the nanoparticle form a two-dimensional vector 
field which, by virtue of the Poincare-Hopf theorem, must 
vanish in some isolated points [15-17 . This topological 
constraint translates into the proverbial statement that 
it is impossible to "comb a hairy ball without creating 
a cowlick". More precisely, one can introduce the no- 
tion of topological charge g of a defect in a vector field 
as the angle (modulo 27r) the vector field rotates in one 
counter-clockwise circuit of any closed contour enclosing 
the defect. Thus a source and a sink have g = 1, while 
a saddle point has q = —1. The Poincare-Hopf theorem 
requires that the sum of the topological charge of all the 
isolated defects (zeros) of a vector field on an oriented 
differentiable manifold M is equal to the Euler charac- 
teristic X of Oil ^ sphere, for which x = 2, any vector 
field must then contain defects, such as a source and a 
sink each with topological charge one (Fig. [2j left). 

The inevitable presence of topological defects on MP- 
NPs has stimulated the imagination of theorists and ex- 
perimentalists who envisioned the possibility of exploit- 
ing functionalized defects as preferential binding sites of 




FIG. 1. Schematic representation of thiol-terminated hydro- 
carbon chains attached to a substrate, (a) Upright configura- 
tion, (b) an optimal tilt and (c) a non-optimal tilt. 



2 




Vertex Au( 100) 

FIG. 2. (Left) Bipolar texture with a source and a sink at 
the north and south pole. (Right) A model of a gold NP with 
octahedral geometry. 



a nanoparticle. In this way one can engineer an entire 
warehouse of "meso- atoms" of controlled valence and di- 
rectionality of bonds [18 . This program has been suc- 
cessfully implemented in the beautiful work of DeVries 
et at. [19] in which gold nanoparticles are coated with 
two species of ligands of different length which then un- 
dergo microphase separation. The defects can then be 
functionalized by chemically attaching sulfur terminated 
chains, such as 11-mercaptoundecanoic acid (MUA). This 
results in a divalent gold nanoparticle with directional 
bonds 180° apart. These in turn can be linked to create 
polymers and free-standing films. The mechanism that 
allows the binding of functional groups to the monolayer 
coating a nanoparticle is referred to as a place-exchange 
reaction and occurs when a thiol-terminated molecule on 
the nanoparticle detaches from the monolayer and is re- 
placed by another molecule from the surrounding solvent. 
As observed in Refs. ^^,20 such a reaction is more likely 
to occur at a spot of the monolayer where a molecule is 
less tightly bound. A defective site is clearly the most 
natural candidate for such a "weak spot", since the tilt 
angle differs here from that favored by vdW interactions. 

The fundamental degree of freedom to describe molec- 
ular tilt on a two-dimensional substrate is the tilt vector 
[21^ corresponding to the projection of the molecular axis 
n on the plane of the substrate: 
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(1) 



with N the substrate normal vector (see Fig. [T]). The 
normalization is chosen in such a way that t = |t | = tan r, 
with r the tilt angle. The interlocking condition between 
chains translates into the requirement t = to = na/d^ 
with a the distance between carbon atoms in a chain and 
d the spacing between neighboring chains j22]. A sim- 
ple description of molecular tilt on coated nanoparticles 
can be achieved by introducing a Ginzburg-Landau (GL) 
model with the following energy functional: 

E = j dA [\KAVitjVH^ + \\(f - tlf] , (2) 

where Vi is the covariant derivative along the surface of 



the substrate. The first term in Eq. ^ favors uniform 
tilt, while the second term implies that molecules have 
a preferential tilt angle tq = arctanto, typically in the 
range 30°- 45° [23H26] . The parameter A is directly re- 
lated to the strength of the vdW interactions between 
carbons. The physically observable states may be de- 
scribed as minimizers of the above GL energy. In the fol- 
lowing we will identify stable equilibrium configurations 
t for the energy (|2| for the case of a spherical substrate 
and an octahedral substrate. The spherical substrate ap- 
proximates the shape of a large NP, while the octahedral 
substrate is appropriate for small NPs where the faceted 
structure of the particle is noticeable (Fig. [2j right). 

Spherical nanoparticles. — On a sphere, the tilt 
vector can be conveniently expressed in a local orthonor- 
mal frame {ei, e2}, namely t = t(cosaei+sinae2). The 
gradient terms can then be expressed as: 

\/^tJ\/H^ = g'^dadjt + t^g'^id^a - A,){dja - A,), (3) 

where g'^^ is the metric tensor and Ai = ei ■ 9^e2 is the 
spin connection, whose curl is the Gaussian curvature of 
the surface: K = e*-^ V^Aj, with e*-^ the Levi-Civita tensor 
[16l [17] . We compute the Euler-Lagrange equations as- 
sociated with Q to find that for a stable equilibrium tilt 
vector field, t and a satisfy the following coupled partial 
differential equations: 
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V- [t\\/a-A)] =0. 



(4a) 



(4b) 



These equations are formally identical to the GL equa- 
tions for superconductors [27j with the t playing the role 
of the absolute value of the wave function, a its phase 
and A the vector potential associated with an exter- 
nally applied magnetic field. Furthermore the quantity 
t^(Va — A) plays the role of the current density and Eq. 
^p) is equivalent to the charge conservation condition. 

As mentioned in the introduction, the simplest config- 
uration formed by tilted molecules on a spherical particle 
consists of two topological defects, a source and a sink, 
placed at the two poles of the sphere (Fig. [2j left). Us- 
ing standard spherical coordinates and (/), it is easy 
to see that, for this bipolar textures, a is constant and 
A = with ^0 = — cos^. Eq. ^p) is then automati- 
cally satisfied provided t is azimuthally symmetric as one 
would expect: 



V • [t^(Va - A)] = -2tA^V^t = 0. 



(5) 



Then, introducing the normalized tilt vector u = t/^o, 
Eq. (|4^) can be written as 

d'^u + cot Odeu + (2/e^ - cot^ 0)u - (2/e^) = 0, (6) 

where = |u| and = 2KA/\tlR^ is a dimensionless 
parameter representing the ratio between the energy cost 
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FIG. 3. (Left) The magnitude of the normahzed tilt vector 
u = |t|/to along a meridian for a bipolar texture for different 
values of e. (Right) The magnitude of the normalized tilt 
vector along the equator as a function of e for the same system. 
The plot shows a transition at e = ^/2. 

due to spatial variations in the tilt of the chains, propor- 
tional to Ka/ (where R is the sphere radius) and that 
of the vdW interactions which favors a uniform tilt and 
is proportional to A^q. 

We display the numerical solution of Eq. (|6| for differ- 
ent values of e in Fig. |3] For small e the magnitude of the 
tilt vector increases rapidly from zero, at the defects, to 
the preferential value t^. On the other hand, for e > a/2, 
the tilted solution is no longer stable and the system un- 
dergoes a transition to a state where all the molecules are 
parallel to the surface normal, and thus t = 0. This sce- 
nario is equivalent to the transition to a superconducting 
phase with the curvature of the substrate playing the role 
of an external magnetic field. Just as the magnetic field 
is given by the curl of the vector potential, the Gaussian 
curvature is given by the curl of the spin connection. The 
bipolar texture is equivalent to the mixed state in which 
the magnetic field penetrates the sample in only a lim- 
ited number of defective spots. Increasing the magnetic 
field eventually destroys the superconducting state and 
the sample is now completely permeable. Analogously, 
at small curvatures (large radii) the chains are uniformly 
tilted with the exception of the two topologically required 
defects at the north and south pole. When the curvature 
is increased above a critical value, the system undergoes 
a transition to a state where t = everywhere or, equiv- 
alent ly, every point is defective. 

Octahedral nanoparticles. — Small nanoparticles 
are far from spherical in shape. In this case we might ex- 
pect the faceted structure to favor the formation of more 
complicated textures that match more effectively the ge- 
ometry of the underlying substrate. Here we consider a 
simple approximation of the actual shape consisting of a 
regular octahedron of edge-length I. The stable equilib- 
ria u satisfy the classic GL-equations on each triangular 
face (i = 1, 2, ••• 8): 

2 

Au+ - (1 -7x^)u = 0, veTi, (7a) 

u = o-(r) , r G ST, , (7b) 

where r the position vector and cr is an appropriate func- 
tion on the boundary dTi of the face T,. Calculating the 
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FIG. 4. Schematic of an octahedral nanoparticle with the 
chains oriented at the vertices of a single triangular face as 
prescribed by the continuous model. 

energy of a given configuration then involves prescrib- 
ing a suitable function a and solving the boundary- value 
problem defined by Eqs. We will in fact consider 

three different boundary-value problems, corresponding 
to different physical scenarios, and prescribe a suitable 
form of cr in each case. 1) The chains rotate smoothly 
across the edge joining two neighboring faces. 2) The 
chains keep their preferential tilt angle along the perime- 
ter of each face. We impose tangential boundary condi- 
tions, which require the tilt vector to be parallel to the 
edges. This necessarily creates discontinuities at the tri- 
angular vertices. 3) An intermediate configuration where 
both tangent and continuous boundary conditions co- 
exist. Both the first and second scenarios have been 
discussed in the literature, where they are sometimes 
referred to as continuous [13] and crystallographic [12] 
models, and using molecular dynamics simulation it has 
been found that they are both possible, depending on 
the length of the chains as well as the temperature of the 
system [14 . 

In the following we will see how these three scenarios 
can be implemented in our framework by a suitable choice 
of a for Eqs. ([t]). The first situation is modeled by 
requiring the chains at the six vertices of the octahedron 
to be parallel to the bisector of the solid angle subtended 
by each vertex (Fig. |4|. Thus looking at the face in the 
positive octant of the lab frame {X, Y, Z}, the chain 
orientations (molecular axes) at the three vertices are 
given by 

ni = X , n2 = Y , ns = Z . 

The normal vector to the face is given by: 

N = (X + Y + Z)/V3. (8) 
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In order to construct the function tr, we introduce a 
smooth interpolation between the three vectors ni, 112 
and Us. The simplest possible choice is a linear interpo- 
lation as shown below: 

1112(5) = (l-s)X + 5Y, 

n23(5) = (l-5)Y + sZ, 

n3i(5) = (1-5)Z + 5X, 

where s G [0, 1] is the normalized arc- length variable 
along each edge, measured counterclockwise, so that 
1112(0) = ni, 1112(1) = n2 and similarly for the vector 
fields n23(5) and n3i(s). The corresponding tilt vector 
can be readily calculated using Eq. ([T]). Expressing this 
in the local x, y frame: 



i(Y-X), y= J|(Z-iX-iY), 



shows that the function a is given by 



(5-l)x+Ji(35-l)y 



sSt- JU3s-2)y 



cri2{s) = 



^23(5) = to 



^31 (s) = 



The magnitude of a is the same on each edge: a = 
Y^2 + 6s{s — l)/to. At the vertices, the tilt angle is thus 
given by Tvertex = arctanA/2 ^ 55°. This depends ex- 
clusively on the geometry of the octahedron and implies 
that any tilted phase whose preferential angle differs from 
55° will be frustrated, in the sense that the preferential 
local order cannot propagate everywhere in the system. 
Other interpolations are also possible, but those do not 
change the overall picture. This configuration is highly 
defective. Because of the natural splay introduced by the 
orientation of the tilt vector at the vertices, each trian- 
gular face has a source at the centroid, a sink at each 
vertex and a saddle at the middle of each edge (Fig. [sj 
left). Hence, there are a total of 26 defects. Calculating 
the total topological charge of this configuration provides 
a realization of the Euler characteristic of convex poly- 
hedra: Qtot = V — E^F = 2, where V, E and F are 
respectively the number of vertices, edges and faces of 
the octahedron and the sign is that of the corresponding 
topological charge. 

For the crystallographic model we impose tangential 
boundary conditions on all three edges. This means that 
on a triangular edge, t has to be parallel to that edge. 
This necessarily creates defects at the vertices: each tri- 
angle has exactly one source, one saddle and one sink lo- 
cated at the vertices. The magnitude t of the tilt vector 




FIG. 5. The tilt vector for the continuous and crystallographic 
models. (Left) Continuous model; the chains smoothly rotate 
across the edges and vertices. This gives rise to a highly 
defective configuration consisting of 8 sources, 6 sinks and 12 
saddles. (Right) Crystallographic model; each triangular face 
contains a source, a saddle and a sink, located at the vertices. 
The resulting configuration has 6 defects. 



will be approximately constant and equal to the prefer- 
ential value to except for a thin boundary layer near each 
vertex where t = 0. A suitable functional form is given 
by t = totanhs(l — s)/e. The function <j can then be 
parametrized as: 



/ N 1 ^(1 ~ ^) - 

^^12 ('5; = tanh x , 



o'2?>{s) = tanh 



.(1-.) 



0-31(5) = tanh ^SL^A (^^^^^f^ ^ 

where e in this case is defined as = 2i^r^/AtQ£^ with ^ 
the edge length of the octahedron. The tilt vector result- 
ing from this model of the crystallographic configuration 
is sketched in Fig. [s] (right). 

The octahedral analog of the bipolar configuration de- 
scribed in the previous section, gives rise to a situation in 
between the continuous and the crystallographic model. 
Using again as a reference the triangle on the positive 
octant of the lab {X, Y, Z} frame, we have that: 



ni 



n2 



ns = N. 



Thus the vector field has a source (sink) at the apex of the 
octahedron, where the chains are parallel to the normal 
vector N, and rotates smoothly across the equator where 
the two square pyramids forming the octahedron join. As 
for the continuous model, we can construct the tilt vector 
along the entire perimeter of a triangular face by linearly 
interpolating between the three vectors ni, n2 and ns. 
This gives us: 

1112(5) = (l-5)X + sY, 
1123(5) = (1-5)Y + 8N, 
1131 (s) = (l-s)N + 5X. 
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FIG. 6. Phase-diagram of the various configurations discussed 
in the text on an octahedral nanoparticle. U represents the 
untilted configuration in which the direction n of the chains is 
everywhere parallel to the normal vector N, thus t = 0. is 
the 6-valent configuration associated with the crystallographic 
model and B is the octahedral analog of the bipolar texture. 



Proceeding as in the previous case, one can construct the 
function cr as shown below: 



0-12(5) = to 



{2s - 1) 



0-31(5) 



to[l + (V3-l)5] 



to[(V3- 1)5 + ^3] 



fx- 



fx- 



With the function cr at hand, we can solve the boundary- 
value problem ([7|) by numerical integration. This was 
performed for different values of e and to using a com- 
mercial finite element software and is summarized in the 
phase-diagram of Fig. [6) where \J represent the untilted 
configuration, 6V\^ the 6-valent configuration associated 
with the crystallographic model and B is the octahedral 
analog of the bipolar texture. As expected, the defec- 
tive configurations are replaced by the untilted one at 
large epsilon, where the energetic cost of spatial varia- 
tions in the chains tilt becomes too expensive. For small 
e both the bipolar and the 6-valent configuration are sta- 
ble depending on the value of the preferential tilt to- The 
smoothly rotating configuration associated with the con- 
tinuous model, on the other hand, is never an energy 
minimizer due to its highly defective structure. Such a 
configuration, however, is likely to occur at finite temper- 
ature where entropy favors the proliferation of defects. 

As noted in the introduction, defects are not only 
mathematically and physically distinguished sites, but 
also they can be chemically functionalized via place- 
exchange reactions. Polymerization chemistry then al- 
lows the functionalized NPs to be linked into chains and 
one can explore polymer physics with polymers based on 



NP monomers. Linear chains of NPs ranging from 3 to 
20 monomers have been observed [19^. Because the lo- 
cal tilt geometry is different for distinct types of defects, 
they are likely to have different chemical reactivity. 6V 
textures with 4 additional defects available for function- 
alizing provide novel possibilities for molecules and bulk 
materials including 2D films, 3D crystal structures and 
branched structures. Observation of any of these higher- 
coordination number structures would provide evidence 
for the texture. Detailed analyses of relative yields of 
various structures could then provide bounds on the rel- 
ative chemical reactivity of different types of defect [28] . 

Conclusion. — In this article we have investigated the 
energetically-preferred configurations adopted by tilted 
chain-like molecules coating a metal nanoparticle by an- 
alyzing a Ginzburg-Landau model for the system. Spher- 
ical nanoparticles have a low temperature ordered phase 
with a bipolar texture characterized by index one defects 
at the poles. Functionalization of such particles together 
with linking gives rise to polymer chains. For small 
nanoparticles the faceted shape of the particle may be 
significant. We extend our analysis to Ginzburg-Landau 
models on octahedra and determine the phase diagram 
as a function of the equilibrium tilt angle to and a dimen- 
sionless parameter e characterizing the relative strength 
of the van der Waals interactions between chains and the 
stiffness modulus Ka controlling spatial inhomogeneities 
in the tilt. We recover an octahedral generalization of the 
bipolar texture but also a qualitatively new texture with 
6 point defects (^V), favored at strong van der Waals 
couplings. Thus shape can dramatically alter the ground 
state and favor extra defects. textures allow for the 
creation of valence 6 NPs after functionalization of the 
defects and consequently new classes of molecules and 
bulk materials when nanoparticles are linked. 
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